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Abstract 



Neutron-proton pairing correlations are studied within the context of two 
solvable models, one based on the algebra SO(5) and the other on the algebra 
SO (8). Boson- mapping techniques are applied to these models and shown to 
provide a convenient methodological tool both for solving such problems and 
for gaining useful insight into general features of pairing. We first focus on the 
SO (5) model, which involves generalized T = 1 pairing. Neither boson mean- 
field methods nor fermion-pair approximations are able to describe in detail 
neutron-proton pairing in this model. The analysis suggests, however, that 
the boson Hamiltonian obtained from a mapping of the fermion Hamiltonian 
contains a pairing force between bosons, pointing to the importance of boson- 
boson (or equivalently four-fermion) correlations with isospin T = and spin 
S = 0. These correlations are investigated by carrying out a second boson 
mapping. Closed forms for the fermion wave functions are given in terms 
of the fermion-pair operators. Similar techniques are applied - albeit in less 
detail - to the SO(8) model, involving a competition between T = 1 and 
T = pairing. Conclusions similar to those of the SO (5) analysis are reached 
regarding the importance of four-particle correlations in systems involving 
neutron-proton pairing. 



21.60.-n, 21.60.Fw,21.60.Jz 



1 



Typeset using REVTgX 



I. INTRODUCTION 



The residual interaction between the nucleons in a nucleus is expected to contain a strong 
neutron-proton pairing component on the basis of isospin-invariance arguments. Practical 
manifestations of neutron-proton pairing have proven elusive, however, in large part because 
most of the nuclei studied to date contain a significant neutron excess. In such nuclei, the 
neutrons and protons near the Fermi surface occupy different valence shells and thus cannot 
effectively exploit the neutron-proton pairing interaction. Moreover, as shown recently [|T],0 , 
even when the active neutrons and protons occupy the same valence shell, they cannot 
effectively build neutron-proton pair correlations except in the very narrow window of iV »s 
Z. 

The development of radioactive beam facilities, now taking place at many laboratories 
worldwide, promises to change the experimental situation dramatically. With these new 
facilities, it should be possible to access all N = Z nuclei up to 100 Sn. In many of the 
heavier N Z proton-rich nuclei, the neutron-proton pairing degree of freedom is expected 
to come into significant play. For this reason, renewed attention is now being devoted to the 
theoretical aspects of neutron-proton pairing. Since full shell-model calculations are feasible 
only for a limited set of nuclei, approximate methods are needed to study in detail such 
collective features. Historically, the method of choice has involved a generalization of the 
usual BCS treatment of pairing between like nucleons |3||f]. Unfortunately, this approach 
does not seem able to provide a correct description of many features of neutron-proton pair 
correlations 0. Isospin projection seems to be a promising avenue to an improved theory 
|U, but up to now it has not been implemented. 

In studying the effects of neutron-proton pairing, simple models can be very useful. 
Several such models, containing a semi-realistic representation of the different modes of 
pairing, have been constructed 0,0]. These models allow for exact solution and simple 
examination of various aspects of approximate treatments. Furthermore, they seem to reflect 
many of the key features of pairing that show up in more realistic shell-model calculations 

In the present paper, we employ the technique of boson mappings || to study pairing 
effects in these models. The basic idea of a boson mapping - as traditionally implemented - 
is to map bi-fermion operators onto boson operators in such a way as to preserve the physics 
of the original fermion problem. In principle, the original fermion problem can be completely 
solved within the boson space. In realistic applications, however, the boson mapping must 
be combined with approximation techniques. 

Boson mappings not only provide a relatively simple methodology for treating the original 
fermion problem, but also can shed light on pairing correlations and on the applicability of 
pair approximations at the fermion level. The bosons are images of fermion pairs. Therefore, 
if the mapping does not result in a simple description in terms of the basis bosons, there 
would correspondingly not be a simple description in terms of fermion pairs. 

Boson-mapping techniques may also be useful in an extended sense, by providing a natu- 
ral methodology for incorporating correlations involving more than just two particles. In the 
neutron-proton pairing models to be discussed, inspection of the boson-mapped Hamilto- 
nian suggests a pair collectivity between the bosons introduced in the mapping. This boson 
pair collectivity can then be treated with a second boson mapping ||, whereby bi-boson 
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operators are mapped onto new bosons representing quartets of the original fermions. The 
original fermion problem can then be rephrased in the language of these new bosons. A 
simple description in terms of these (quartet) bosons would confirm the importance of the 
associated four-fermion correlated structures. 

The paper is organized as follows. We first consider the SO (5) model of monopole- 
isovector pairing. In Section U we briefly discuss the model and its boson mapping and 
then examine in detail several variants of boson mean-field approximations. The applicability 
of fermion-pair approximations is then studied in Section [HI]. In Section [IV], a second boson 
mapping of the SO(5) model is performed and the role of bi-boson (four-fermion) structures 
is investigated. We then turn to the SO (8) model with both isoscalar and isovector degrees 
of freedom in Section 0. The same issues are addressed as for the SO(5) model, but without 
as much detail. Finally, Section [VI| summarizes the key conclusions of the work and spells 
out some issues for future consideration. 



II. SO(5) MODEL - BOSON MEAN-FIELD METHODS 

A. The SO (5) model and its boson realization 

The SO (5) model is perhaps the simplest tool to meaningfully investigate neutron-proton 
pairing. In this model, a system of N nucleons occupying a set of degenerate single-particle 
orbits with total degeneracy 4f2 = 2J2{2j + 1) interact via an isovector pairing interaction, 

H = grf-S, (1) 

where 

sJ^EjV'^IE , (2) 

Z 3 

and 

s v = {-y- u s_ P . (3) 

The Hamiltonian ([I]) is invariant under the group of SO (5) tranformations generated 
by the three pair creation operators Si, the three conjugate pair annihilation operators S u , 
and the three components of the isospin operator T. This makes the analysis of the model 
extremely simple and it is for this reason that it has been used recently in several studies of 
relevance to neutron-proton pairing 



A Dyson boson realization of the SO (5) algebra has been constructed in Ref. [|1(|. For 
the isovector pairing model, the boson space is constructed in terms of a scalar (L = S = 
J = 0) isovector (T = 1) boson s„, and the mapping takes the form [[| 

Sy ^ S v , 

% - V2[s^sZ . (4) 
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Here, N = —s^ ■ s is the boson number operator. In these and all subsequent expressions, 
we use the standard notation for the scalar product, 



* f • * f = E(-r4s f - 



B. Approximate boson mean-field methods 



Boson mappings provide an alternative technique for solving a fermion problem. Diag- 
onalizing the mapped Hamiltonian in the ideal boson space (where possible) would yield 
all of the eigenvalues and the boson images of all eigenvectors of the original fermion prob- 
lem. Note, however, that in the boson space spurious states may occur. These are boson 
states with no counterparts in the original fermion space that arise as a pure artifact of the 
mapping. In physically realistic problems, where it is impossible to diagonalize the mapped 
Hamiltonian in the full boson space, these spurious states cannot be readily removed from 
the problem. In the SO (5) model, however, this is not the case. Here it is possible to exactly 
diagonalize the mapped Hamiltonian. Furthermore, it can be shown that the spurious states 
only arise for N > 2Q and that they all have isospin T > 2Q — N/2. Thus, by focusing 
our analysis on systems with iV < 2Q, we can be sure that spurious boson states do not 
contaminate the physics of interest. 

In the present study, we also treat the dynamics of the boson problem using approxi- 
mate methods. By comparing with the exact results, we can assess the usefulness of these 
approximate methods in capturing the dominant collective dynamics in the presence of pair 
correlations. 

We begin by considering boson mean-field techniques. The starting point is to introduce 
a mean-field boson 7 as a linear combination of the basis bosons. The creation operator for 
the mean-field boson can thus be expressed as 



From this mean-field boson, several approximate variational states can be considered. 
Minimization of the energy of these variational states subject to the relevant constraints 
defines possible boson mean-field approximation procedures, each of which is the natural 
analog of a fermion variational procedure. We will denote each boson variational procedure 
by the corresponding standard fermion terminology. 

In the boson analog of the BCS method, for example, the variational state is 
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t _ 



7i«i + 7-iS-i + 7oSo • 



(5) 



|BSC) cxexp(?77 t )|0) , 



(6) 



with constraints on the number of bosons (one-half the number of fermions) 



(BCS|iV/2|BCS) = Af 



and on the ^-component T z of the isospin 



(BCS|T |BCS) = T z 
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The boson analog of number-projected BCS is the Hartree-Bose procedure, whereby the 
variational state is of the form of a boson condensate^ 



|BSC,A0 oc7 tA/ ]0) (7) 

with the single constraint 

(BCS,A/1T |BCS,A0 = T z . 

The next level of mean-field approximation we consider involves the number- T z projected 
BCS variational state, 

\B$C,MT z )ocV T ^\0) . (8) 

Here, Vt z projects the boson condensate onto a state with definite T z . 

Finally, if we were to apply, in addition to number and T z projection, full isospin T 
projection, the exact SO (5) state would be realized (as long as the probe function has a 
nonzero overlap with the exact state). 

C. Energies 

We now compare the energies arising at the various levels of boson mean-field approx- 
imation with the exact ground-state energies of the SO(5) Hamiltonian ([!]). We consider 
systems with an even number of nucleons only 

The exact eigenenergies for a system with M nucleon pairs are given by 

E = -g[{N-\v s ){n+\-\N-\v s )-\T{T + l)\ , (9) 

where v s is the singlet-pairing seniority 

We are especially interested in the ground state of the system, which is realized for 
v s = 0. The ground-state energy is then given by 

^exact = - 9 mn+^-^) 

-l Tz (T z + l)-8(T z + l)} , (10) 

where 5 = for even-even (T = T z ) systems and 1 for odd-odd (T = T z + 1) systems. 

Applying the Dyson boson mapping (^) to the SO (5) Hamiltonian (|I]) leads to the boson 
Hamiltonian 

H B = -g[Ar(Q + l-AT) + ^ ■ s^s-s] . (11) 



Throughout the present paper, we understand the projection as performed before variation. 
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When this Hamiltonian is used in the variational approximations described in the preceding 
subsection, only two kinds of solutions can occur. The first, which we denote as A, has 
7o = 0; the second, denoted by B, has 71 = 7_x = 0, 70 = 1.0 

In the BCS case, with trial state @, the solutions A and B are degenerate for T z = 0, 
whereas only solution A applies for T z > 0. Note that this is precisely what was found in 
the generalized BCS treatment of Ref. Q. The variational energy for these solutions is 

£(BCS) = £ exact + g[^Af - ^T z - 5(T Z + 1)] . 

Adding number projection via the trial state ([?]) leaves many properties of the solution(s) 
unchanged. The solutions A and B are still degenerate for T z = 0. And, for T z > 0, only the 
solution A applies. There is a change in the variational energy, however, which now becomes 

£(BCS,A0 

= £ exact + gM\l - + ?*)] - gS(T z + 1) . (12) 

Note that for T z — > Af, the approximate energy at this level of approximation approaches 
the exact result, as expected for a degenerate-orbit pairing model with one type of nucleon. 

When imposing T z projection as well, some new features appear. For even-even nuclei, 
for example, T z projection removes the degeneracy of the solutions A and B, the solution A 
giving the lower energy. Thus, for even-even nuclei, the ground state solution is of the form 
A with variational energy 

E(BCS,AfT z ) = E exact + g^Af(l - ^) . (13) 

Here, too, for T z — > Af, the approximate result approaches the exact result, as clearly it 
should. 

For odd-odd nuclei, on the other hand, the solution A disappears. Clearly, an odd-odd 
system must contain at least one T z = (sq) boson, and such components are not present 
in solution A. 

In fact, problems with the solution A in odd-odd nuclei already show up at the level of 
number-projected BCS approximation, though not as transparently. For an odd-odd nucleus 
with T z = Af — 1, the number-projected variational energy (see Eq. 12) is lower than the 
exact energyi 

For the odd-odd T z = nucleus, the number-projected solution is thus of the form B and 
already has good T z . Furthermore, the variational energy at this level of approximation is 

E(BCS,AfT z = 0) = £ cxact + g(Af - 1) . 



2 This notation is in correspondence with Ref. 0]. 

3 The same effect is seen in the fermion generalized BCS calculations of Ref. 0]. 
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Tboson mean field 



ferniion-pair approx. 




FIG. 1. The exact and BCS S0(5) ground-state energies E in units of g, two-particle transfer 
strengths S nn and S np , the Fermi strengths Sp" , and overlaps of exact and BCS wave functions 
are shown as a function of T z for 0, = 10 and N = 10. The results of the boson mean-field 
approximations and the fermion-pair approximations are displayed in the left and right parts of 
the figure, respectively. The exact results are denoted by circles, the BCS results by squares, the 
number projected BCS results by triangles pointing up, and the number-T 2 projected BCS results 
by triangles pointing down. 
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For the odd-odd T z > nucleus, there is no number-T 2 projected variational solution 
from the above- discussed class of trial functions. 

A comparison of the variational energies obtained at the various levels of approximation 
discussed above with the exact ground-state energies is given in Fig. la. Two points 
are immediately evident from the figure. First, as expected, the successive restoration of 
symmetries improves the agreement with the exact results. Second, the variational energies, 
independent of the specific mean-field approximation, agree quite well with the exact values. 
This can be readily understood from the structure of the boson Hamiltonian ([11]). In this 
Hamiltonian, the first term only depends on J\, and its contribution to the variational 
energy for states with a given number of bosons is not dependent on the form of the wave 
function. Moreover, the largest contribution to the variational energy from the second term 
— ~gs^ ■ s^s ■ s in ( |TTj) is —^gAf 2 for all the above- discussed solutions. As a consequence, 
the ground-state energy is not a particularly sensitive observable for assessing the quality of 
the various approximate methods, at least for this model. 



D. Two-particle transfer strengths 

To assess the quality of the different variational solutions, we must consider other observ- 
ables as well. In the SO (5) model, there are not many such possibilities. The pair operators 
multipled by 1/ y/Ti give the normalized two-nucleon transfer strengths of the respective pair 
PH . The mean values of the product of pair operators may then be connected with the 



summed strengths for two-nucleon pick-up from a given state i [12 



S np = (Ar,i\^SlS Q \Af,i) 

= J2\<Af-l,f\-Ls \Af,i)\ 2 , (14) 



/ 

and similarly for S nn and 5* pp .i Since two-neutron-transfer reactions provide valuable infor- 



mation on nn pairing correlations [13], it is expected that np transfer will analogously be 
sensitive to the neutron-proton-pairing mode. The sum of the three total transfer strengths 
is simply related to the energy of the spin-isospin conserving Hamiltonian ([!]). 

Exact values of the total transfer strengths for each mode can be readily deduced from 
formulae given in Refs. [ttlH. For the approximate boson mean-field methods, we find 



BCS: 



1 \f + T 

^(solution A) = -{M + T Z ){1 



2 V /v 2Q 

^(solution B) = JV(1 - ^r) . 



4 These mean values are related to the pairing gaps of the respective pairing modes ||. Their 
usefulness as a rough measure of the number of the respective pairs has also been extensively 
discussed ggg. 
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• Number-projected BCS: 

S nn (solution A) = -(N + T z )(l - 

S np (solution B) = JV(1 - 

• Number-T 2 projected BCS: 

5 nn (solution A) = ^{AT + T z )(l 

S np (solution B) = Af(l - 

• All methods: 

S nn (solution B) = , 

S np (solutionA) = . 

The exact and aproximate results for the pair-transfer strengths are compared in Fig. 
lb-c. The exact and approximate values differ considerably, much more so than the energies. 
For example, for T = states, all three exact total transfer strengths are equal JO. The 
approximate methods, however, give a zero transfer strength for the mode that is not present 
in the approximate wave function. Thus, even though the sum of all transfer strengths - as 
mirrored in the energy - is quite close to the exact value, its composition from the various 
terms of the Hamiltonian may be completely wrong. 

E. Fermi strengths 

Two other physically-interesting observables in the SO (5) model are the Fermi total 
strengths 

S£ =< T~T + > , 

S F =< T + T- > . 
These two quantities are related by the Ikeda sum rule 

S-p Sp = 2iT z . 
The exact expression for the Sp strength is 



M - lM + T Z 
AT 2VL 



A/--1. 
2tt ' 



Af + T z -2^ 
2VL 



AT- 1 
2VL 



) ■ 
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Spt(exact) = 25{T Z + 1) . 

The results for the two types of approximate solutions, A and B, are independent of the 
specific boson mean-field method. They are given by 

(solution A) = M - T z , 



^(solution B) = 27V . 

A comparison of exact and approximate results is given in Fig. Id. We see that the 
approximate results may differ substantially from the exact ones, especially for small values 
of T z . This again suggests that even if the approximate wave functions provide reasonable 
values for energies their quality may be quite poor for other more-sensitive observables. 



F. Overlaps 



Another measure of the quality of approximate methods is the overlap of the approxi- 
mate wave functions with the exact ones. Such a measure is not particularly useful when 
considering the BCS and number-projected BCS approximations, however, where the ap- 
proximate wave functions are averaged over nuclei with different number of nucleon pairs 
and/or different T z . Therefore, we calculate the overlaps only for the number-T^ projected 
BCS method. The results are illustrated in Fig. le. 

We see that for small values of T z , the approximate solution has very little overlap with 
the exact solution. It is only for T z — > M that the overlap approaches one. In this limit, 
the solution A (a pure nn condensate) is the exact ground state of the model system. In the 
limit of symmetric nuclei, however, the approximate ground state wave function is very bad 
and a more sophisticated procedure is necessary in order to get an acceptable description of 
the system. 



III. SO (5) MODEL - FERMION-PAIR APPROXIMATIONS 



In this section, we discuss the solution of the SO (5) model in the original fermion space. 
While the exact solutions and the BCS solutions have been discussed before, the effect of 
number and T z projection in the BCS approach has to our knowledge never been studied. 

As noted in Subsect. [11 B|, each of the boson mean- field approximations that we consider 



is the analogue of a fermion-pair approximation. The variational wave functions associated 
with these fermion-pair approximations can be readily obtained from (|J), (^]), and ([5]) by 
replacing the collective mean-field boson creation operator 7"^ with a collective fermion-pair 
creation operator 

r+ = r 1 s , 1 t + r_ 1 s , i 1 + r sj (15) 

and the boson vacuum |0) with the fermion vacuum |0). 

Despite the obvious similarities between the boson mean-field and the corresponding 
fermion-pair approximations, they are not identical. Pauli effects are accomodated very 
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differently in the two approaches and this can lead to differences in the results. It is only in 
the limit Q — > oo that the corresponding bosonic and fermionic results agree. Nevertheless, 
we expect the main features found in the boson mean-field analysis to be present in the 
fermion-pair methods as well. 

Such a conjecture is readily confirmed for the SO (5) model under investigation. In the 
fermion-pair approximations, for example, two types of solution likewise occur. There is a 
solution A corresponding to the nn — pp phase (To = 0) and a solution B corresponding to 
the np — pii phase (Ti = r_i = 0). Furthermore, as in the boson mean-field treatment, (i) 
T z projection removes the degeneracy of the solutions A and B, and (ii) both solutions are 
unphysical for odd-odd nuclei, except in the case of T z = where the solution B is relevant. 

Simple analytic expressions can be obtained in the fermion-pair approximations for many 
of the quantities discussed in Sect. [TI|. In the fermion BCS, for example, the following simple 
expressions obtain: 

£(BCS) 

= ^exact + gW{~ ~T Z (1 + ^L) - 5(T Z + 1)] , 

1 J\f + T J\f + T 

^(solution A) = UN + T z )(l - 



2 V Jy 2Vt 2Vt 2 



M 2 

S nn (solution B) " 



4Q 2 ' 



J\f 2 - T 2 

S np (solution A) = — ^ 

AO? 



S np (solution B) = jV(l - ^- + A 



^(solution A) =Af -T z 



2tt 40 2/ 

M 2 - T 2 
2VL 



^(solution B) = 2Af(l - ^) . 

For the number and number-T^ projected BCS, we have not obtained closed expressions 
for all of the quantities of interest, even though it is likely that they too can be derived. 
For those cases, the results we present have been obtained numerically. The one quantity 
for which we have obtained an analytic expression is the energy of the number-T z projected 
state, 

a r Af 2 — T 2 

E{BCSMT Z ) = £ cxact + |/V(1 - ^ - ) • (16) 
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In the right half of FigfJ, the results of the approximate fermion-pair methods are shown 
next to their corresponding boson mean-field results. The fermionic results are in general 
closer to the exact ones than their bosonic counterparts. Nevertheless, the main deficiencies 
found in the bosonic analysis persist when working directly in the fermion space. Most 
notable are the large disrepancies in the Fermi and two-particle transfer strengths and in 
the overlaps between the approximate and exact wave functions. 

Closer inspection of the results indicates that even for the relatively well-reproduced 
ground-state energies the approximate methods do not capture some important details. An 
example is the double binding energy difference [J14||, 



SV np (N, Z) = ±{[B(N, Z) - B(N - 2, Z)\ 

— [B(N, Z - 2) - B(N - 2, Z - 2)]} . 

For even-even nuclei in this model, the exact expression for this quantity is 

5V np (N,Z) = ~g ,{orN = Z 
= , otherwise . 

The jump of 5V np (N, Z) at the N = Z line, which is a persistent feature of experimental 



data [jlJl, may be used to isolate the Wigner term in the binding energy ||16|| .13 

When 5V np (N, Z) is calculated using the fermion number-T z projected BCS method 
- the approximation that yields the best reproduction of the exact energies - one obtains 



5V ap (N,Z) = --^g . 



irrespective of whether N = Z or not. Clearly, the physics of the jump at N = Z is rather 
subtle and its correct description needs more sophisticated approximate methods than BCS 
or any of its variants. 

Apparently, the standard fermion-pair approximations do not allow for the coexistence 
of like-particle and neutron-proton pairs in the case of an isospin- conserving Hamiltonian. 
This results in a very small value of the two-nucleon transfer strength for the mode not 
present in the wave function. As shown recently ||, all three pairing modes can coexist, 
however, when the isospin symmetry of the Hamiltonian is violated. This led to a proposal 
to introduce an isospin-breaking Hamiltonian, chosen to reproduce at BCS level the various 
pairing gaps in symmetric nuclei, and to use it in BCS treatments of all nuclei. 

We question such a prescription. One can rephrase the prescription of Ref. [||] and 
look for the parameters T u in the general pair creation operator fll5|) that in T = nuclei 
give appropriate two-nucleon transfer strengths (or pairing gaps) for all three modes. The 
condition within the SO (5) model is that all three strengths should be equal. Working in 
the number-T z projected BCS approximation and focusing on the case Q — 10, H — 10, we 



Of course, the actual Wigner term may not have its main contribution coming from isovector 
pairing, as it does in the SO (5) model. 
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find that this equality can be achieved by choosing = r 2 1 = 0.2593, and Tq = 0.4815. 
When we then calculate the ground-state energy and the Fermi strength with this same 
choice of structure coefficients (for the same T = nucleus), we obtain the results -10.71 
and 108.6, respectively. These are to be compared with the exact results of -65.0 and 
0, respectively, calculated using the isosp in- invariant Hamiltonian. Clearly, the proposed 
prescription cannot describe in a unified way the various observables of interest. 

We can understand these results as follows. The fermion-pair state just described, namely 
the one that yields equal pair-transfer strengths for all three modes, is in fact very close to 
the exact state with maximal isospin T max = Af, 

\BCS : ATT max T,) cc V Tz {\s\ + \sU + 4=Sjy» . (17) 

This suggests that the isospin-breaking procedure proposed in Ref. || produces high-isospin 
admixtures in the ground state that are too large. And this makes questionable its usefulness 
in describing the ground state of an isospin-conserving Hamiltonian when T z rs 0. 

It is useful to summarize here the principal findings up to this point in the analysis. It 
has been argued many times that the (generalized) BCS method cannot properly describe 
neutron-proton pairing in the SO (5) model ]|]. The present results show that this conclu- 
sion does not change when number and T z projection are switched on (with variation after 
projection). The deficiencies of these standard fermionic methods for treating pairing are 
clearly seen in the analysis. And, furthermore, they can be alternatively seen within the 
context of analogous mean-field boson methods applied following a boson mapping of the 
model. 

IV. SO (5) MODEL - BEYOND FERMION-PAIR CORRELATIONS 

Since fermion-pair approximations and the corresponding boson mean-field methods 
show deficiencies when applied to the ground state of the isovector neutron-proton pair- 
ing SO (5) model, we need to consider more sophisticated procedures. 

A. BCS for boson Hamiltonian 

The boson-mapped SO (5) Hamiltonian ( |TT| ) contains, in addition to the dominant linear 
term, an attractive boson-boson pairing interaction. Two bosons only interact when in 
a J = T = scalar-isoscalar state. A natural first approach to consider for a system 
dominated by boson pairing is boson BCS approximation fT7j| . 

In this approach, a variational boson wave function of the form! 

|$) oc exp(r/4 + Cs f • s f )|0) (18) 
is considered. A generalized Bogolyubov transformation, 



'We limit our discussion here to even-even systems. 
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a l = us l - vs v - x5i u , (19) 
is then introduced, with the constraint 

2 2 1 
U — V = 1 . 

The state flT8| ) is the vacuum for quasiboson operators flI9| ) if the relations 

2 M ( + u = , 
u?7 — x = 

are satisfied.! Constraints on the average number of pairs N and the average value of the 
isospin T z give 

2 M-T z 
v — 



2T Z + 3 
x 2 = T, . 

The energy of the state ( |18D is then easily calculated. We present explicitly the result for 
T z = only, for which the energy is 

£(BCS) = £ exact + Af(l + |A0 • (20) 

An illustrative comparison of the exact and approximate boson BCS energies is given in 
the upper part of Fig. |2|. We see that the boson BCS theory cannot explain the energies 
satisfactorily. 

Its failure can be traced to a very large dispersion in the pair number contained in the 
wave function fllTf). For the model under discussion, the dispersion is given by 



m? = ( ,^|.) - * = * ~ T >)^ ; + ^ + «) + T , . (21 ) 

The results are displayed in the lower part of Fig. |2|. As is clearly evident, the dispersion is 
quite large, especially for small values of T z . 



In fact, the wave function (y) of Subsect. II B can be considered in an analogous way with v = 0. 
We denote the method in the present subsection as the boson BCS to distinguish it from the BCS 
in the mean-field boson approximation. In the latter, the BCS bra state is a boson image of 
the fermion BCS state. 
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FIG. 2. The exact and boson BCS SO(5) ground-state energies in units of g are shown in the 
upper part of the figure as a function of T z for f2 = 10 and N = 10. The pair-number dispersion 
in the boson BCS state is given in the lower part. 

Some understanding of this result can be obtained by focussing on the case of T z =0. 
Then, Eq . (|2l|) is a special case of a formula for the dispersion associated with the BCS wave 
function ( |T8| ) for rj = and for pairing in a single level of degeneracy 2Q: 

(22) 



(AAO 2 = M 2 /Q + 27V . 
In the case under discussion of an isovector s boson, 2Q = 3 and 



(AM) 2 = -Af(Af + 3) . 



(23) 



Because of the inverse dependence on f2 of the dominant quadratic term in (22]), the disper- 
sion is quite large at T 2 =0.i 

Thus, due primarily to the small dimensionality of the s-boson space (d = 3), boson BCS 
is not an acceptable approach for treating boson-boson pairing correlations in this model. 
It should be noted, however, that when number and T z projection are applied to QT%D the 
exact ground-state results are obtained. 



Note that in boson BCS theory, as contrasted to fermion BCS theory, the linear and quadratic 
contributions add coherently, giving a further reason for the large dispersion. 
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B. A second boson mapping 

As we saw in Section 2, a boson mapping followed by a mean-field treatment of the 
resulting Hamiltonian is an alternative to a BCS treatment of the original problem. Thus, 
we will now consider the possibility of implementing a second mapping of the problem, from 
the system of s v bosons (representing fermion pairs) to a new system of bosons that represent 
fermion quartets. We will first develop and apply the ideas to systems with T = and then 
subsequently (for reasons to be clarified later) to systems with T ^ 0. 

1. T = case 

In the boson- mapped SO (5) Hamiltonian ( |TT|) only the operators 

■ , s-s , s f ■ s (24) 

appear. This set of operators closes under commutation, 

-2s ] ■ s ] , 
2s -s , 

-4^-5 + 6 , (25) 

generating the algebra 0(2,1). 

In the second boson mapping, the boson space built up from the three s v bosons is 
mapped onto a new space built up in terms of a single J = T = boson, which we denote 
t. Whereas the s u bosons represented pairs of the original fermions, these new bosons 
represent pairs of s bosons, or equivalently quartets of the original fermions. It should be 
clear at this point that the mapping just outlined can only give information about states 
with total T = 0. 

In more detail, the mapping is constructed such that the commutation algebra of bi- 
boson operators ( |2"4"D in the s-boson space are preserved in the t-boson space. The Dyson 
realization of the mapping for the 0(2,1) algebra is 

s f • s f -> 6t f + Athh , 
s ■ s — * t , 

s f ■ s -> -2tH . (26) 

Here, v and t denote the creation and annihilation operators of the boson t. 
We obtain after the second boson mapping a U(l) Hamiltonian 

H BB = -2g(n + | - th)th . (27) 

The eigenfunctions and eigenvalues of this Hamiltonian can be simply obtained. The 
eigenfunctions are condensates of t bosons, |t n )) i and the corresponding eigenenergies are 



9 We use a double bracket )) to denote a state in the £-boson space. 



W ■ s, s' 



s.s ■ s\ 



s ■ s 
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E = -2g{(Q + ^-n)n . (28) 

These eigenenergies are exactly equal to the ground-state energies (|T0|) of the original fermion 
SO (5) Hamiltonian (JJ) if the obvious relation n = |/V is invoked. 

There are several points that should be made before proceeding to the T ^ case: 

• The exact eigenstate \t n )) is in the form of a Hartree-Bose approximate solution for 
the Hamiltonian i^BB- Equivalently, following the terminology of Sect. (HO, we can 
think of it as representing the number-projected BCS approximation of H-q. Thus, 
we have confirmed that BCS approximation of the mapped Hamiltonian H-q is an 
acceptable way to describe the full correlation structure of the problem, as long as 
number projection is included. 

• Clearly, the solution |t n )) involves four-fermion correlations since the boson t itself 
represents a correlated quartet. Thus, we have confirmed the essential role played 
by four-body correlations in systems involving both like-particle (nn and pp) pairing 
correlations and np pairing correlations. 

• Lastly, we have demonstrated, through the example studied here, the possible use- 
fulness of iterative boson mappings in accomodating many-particle correlation struc- 
tures. Considering the difficulty in building quasiparticle methods that can treat 
many-particle (e.g., quartet) correlations directly at the original fermion level, we feel 
that this could be an important conclusion. 

Of course, we should bear in mind that all of the above points have been demonstrated so 
far for T = systems only. 



2. T ^ case 

The boson space introduced in the previous subsection is constructed entirely in terms 
of the J = T = t boson. Thus, as noted earlier, the Hamiltonian fl2"7D derived using the 
mapping ( ffif) to that space can only provide information on even-TV J = T = states. 

To make possible a more complete analysis, we have to extend the set of operators (p4[) 
considered for the second boson mapping. If we want to also include nuclei with M = odd 



and/or T z 0, we have to add to the set of operators (|24j) the creation and annihilation 
operators of the s boson, 

s f , s . (29) 
We should thus consider, in addition to (|25[), the commutation relations 



[s v , s T ■ s T ] = -2s u , 



t . c tl — _9et 
, 

-s v , (30) 



and their hermitian conjugates. 
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The algebra of the operators (p4| ) and ( p29|) have their Dyson boson realization on the 
space formed by the scalar-isoscalar boson t (representing,as before, fermion quartets) and 
the scalar-isovector boson a v (representing fermion pairs) .1121 

The Dyson realization of this extended algebra is 

s f • s f -> 6t f + tthh + a ] ■ a ] - tt ] a ] ■ a , 
s ■ s — > t , 

s f • s ->• -2t f t + a ] ■ a , 
4 -»■ <7+ - 2tt^ , 

s u — > a v . (31) 

The ideal space formed by the t and a bosons is larger than the original s-boson space. 
Nonphysical (spurious) states are introduced by the second boson mapping, in addition to 
those already introduced by the first boson mapping. Full diagonalization of the mapped 
Hamiltonian in the ideal space separates physical states from spurious states, however, pro- 
viding eigenvalues and wave function images of the original problem. 



Applying the second boson mapping ([H]) to the boson Hamiltonian (|TT|), we obtain 



H BB = -g[Af(Q + 1 - A)") + 3th + 2thht 

-2tfto* ■ a + -a f ■ ah] , (32) 

where we have used the notation 

fif = 2th -a ] - or 

for the total number of nucleon pairs in the system. Note that the Hamiltonian (|32|) , when 
restricted to states built from t bosons, correctly reduces to the T = Hamiltonian (|27|). 
Though the Hamiltonian (32) is nonhermitian (a general feature of Dyson mappings) 



with non-zero nondiagonal elements, we can nevertheless easily determine its eigenenergies. 
In the basis with good isospin T, 



the Hamiltonian matrix for (^) has all zero elements below the main diagonal. Therefore, 



the diagonal matrix elements coincide with its eigenvalues, which we find to be 

-g[M(n + 2 - iAO - \n a {n a + 1)] . 



°This procedure is analogous to the familiar boson- fermion mappings. In those mappings, intro- 
duced to simultaneously describe systems with an even number and an odd number of fermions, the 
set of bi-fermion and fermion operators is mapped onto a space formed by boson and quasi-fermion 
operators. 
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The lowest eigenvalue, with n a = T, corresponds to the physical state. O The other eigen- 
states are nonphysical. The energies of the physical solutions are in precise agreement with 
the exact energies (0) of the v s = states, as they should be. 

The corresponding left eigenvector for the physical n a = T state is readily found to be 

((t^- T ^;a T T\ . (33) 

Note further that the exact left eigenstate is a product of a t-boson condensate and a a-boson 
condensate, precisely the form of a coupled Hartree-Bose solution. 

From the above remarks, we see that all of the conclusions given at the end of the 
preceding subsection (for T = states) carry over to the more general case. Namely, 
(i) number-projected boson BCS approximation following the first boson mapping is an 
appropriate method for incorporating the correlations contained in the SO(5) model, (m) 
incorporating four-particle correlations within the context of alpha-like clusters is the key to 
describing the structure of the system, and (Hi) iterative boson mappings are an attractive 
means of accomodating the various correlations contained in the model. 



C. Wave function in terms of fermion-pair operators 

In the previous subsection, we obtained an analytic form for the left eigenstates of the 
SO (5) model in the t — a space of the second boson mapping. For this particular problem, 
we can invert the procedure to obtain the corresponding exact fermion wave functions. 

The key to the procedure is to focus on the left eigenstates ( |3~3"D of the second boson 
Hamiltonian. As can be seen from the mapping equations (|4], ^ and |31]), there is a simple 
chain of relations that take us from the left fermion eigenstates to the left second boson 
eigenstates, namely 

S ■ S — > s ■ s — > t , 

S u ^s v ^o v . (34) 

We can therefore write the left fermion eigenstate that leads to (|3^) asll 

(ATT\ oc (TT\(S ■ S)^-^ . 

Since at the fermion level the right and left eigenstates are conjugates to one another, 
we can now write the right eigenstate as 

\NT) oc (§■ S)^ T) \TT) . 

Finally, making use of the relation (|ITD , we obtain for the exact SO (5) state vectors in 
terms of fermion pairs 



11 For J\f > fl, the condition T < 2fl — J\f must also be taken into account. 

12 Note that the normalization of the fermion state cannot be simply deduced from the boson 
state, so that only a proportionality relation is given. 
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WTT Z ) 

oc P Tm (& ■ S^- T \\S\ + \SU + ^4f\0) . (35) 

A special case of this relation for the ground states of even-even nuclei (with T = T z ) and 
odd-odd nuclei (with T = T z + 1) was given in Ref. 0. Note that the exact wave functions 
involve a condensate of ■ J=0 T=0 quartets and an isospin-stretched condensate of 
J=0 T=l pairs projected to good T z . 

The above considerations reconfirm what was demonstrated in the previous subsections, 
namely that four-particle correlations are an essential ingredient for a correct description of 
the eigenstates of the neutron-proton pairing SO (5) model. 

V. SO(8) MODEL 

A. The model and its boson realization 

A second algebraic model that has been used extensively RPlfl to study neutron-proton 
pairing correlations is one based on the algebra SO (8). As in the SO (5) model, neutrons 
and protons move in a set of degenerate single-particle orbits of total degeneracy 4Q. In 
this model, however, they interact via a sum of an scalar-isovector pairing interaction and a 
vector-isoscalar pairing interaction, with the Hamiltonian taking the form0 

ff= ?(l±f) 5 t.s + s(l^i)p».p +9ph ^.^. (36) 

Here, the operator Si creates the same L=0 S=0 J=0 T—l isovector pair as in the SO(5) 
model, the operator creates an L=0 S—l J—l T=0 isoscalar pair, and JF^ is the Gamow- 
Teller operator. 

Note that the relative importance of isoscalar and isovector pairing in the Hamiltonian 
(j36|) is governed by a single parameter x, which varies from —1 (pure isoscalar pairing) to 
+ 1 (pure isovector pairing). The last term in the Hamiltonian, a particle- hole force in the 
T = 1 S = 1 channel, is included to bring the Hamiltonian into closer contact with more 
realistic nuclear Hamiltonians, without destroying the simplicity of the model. 

The Hamiltonian (|36D is invariant under the group of SO (8) transformations generated 
by the three isovector pair creation operators Si, their three conjugate annihilation oper- 
ators S u , the three isoscalar pair creation operators P^, their three conjugate annihilation 
operators P M , the three components of the isospin operator T v , the three components of the 
spin operator iS M , and the nine components of the Gamow- Teller operator T^. 

The Dyson boson realization of the SO (8) algebra is constructed by mapping its bi- 
fermion operators onto bosonic operators formed from the creation operators si of a scalar- 
isovector boson s and pj, of a vector-isoscalar boson p, and their conjugate annihilation 
operators s u = (— ) 1 ~ u s_ u and = (— ) 1_M P- At 
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We only consider the isospin-conserving version of the SO(8) model in this work. 
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TV" = —(s^-s+p^-p) 



(37) 



is the total boson number operator. 

As in the SO (5) model, spurious states are introduced by this mapping. Here too they 
arise for N > 2Q and can be readily identified when a boson basis with good spin and 
isospin is employed. 



B. Dynamical symmetries of the model 

For certain values of its parameters, the SO (8) Hamiltonian (^) exhibits dynamical 
symmetries. In such cases, the eigenvalues can be obtained analytically. For all other 
parameters, analytic expressions for the eigenvalues cannot be obtained, and numerical 
diagonalization is required. 

The SO T (5) dynamical symmetry limit of the model is realized when x = 1 and g p h = 0. 
In this case, the Hamiltonian reduces to (p]) and its ground-state energy is given analytically 
by©- 

The SO s (5) limit is realized for the parameters x = — 1 and g p h = 0. In this case, the 
ground-state energy is given by 

E = -g[(M -T z )(tt + ~- ±M -±T Z ) - 8} . (38) 

The third dynamical symmetry arises when x = 0, namely when there are equal amounts 
of isovector and isoscalar pairing. In this case, the Hamiltonian fl3"ED has an SU(4) dynamical 
symmetry and its exact eigenvalues can be written as 

E = --g[2H{n + 3) - M 2 - A(A + 4)] 

+g ph [X(X + A)-S(S + l)-T(T+l)] , (39) 

with A the usual SU(4) label. 

The ground-state solution in this symmetry limit arises when A = T z + 5, and the 
corresponding energy is 

E = --g[2Af{n + 3) -Af 2 - T Z (T Z + 4) - 5(2T Z + 5)] 

+3<? ph (T z + 5) . (40) 
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C. The first boson-mapped Hamiltonian 

Mapping the fermion Hamiltonian flSop onto the space of s and p bosons leads to the 
boson Hamiltonian 

h b = -£0d_^[(n + 1 - a>V . g + I( p t . p t _ a t . s t)s . §)] 
- ^ 2 [(» + 1 - • P + ■ s f - p ] ■ p ] )p ■ p)\ 

+g p h\p ■ P^s ■ s + • s'p ■ p + 2s* ■ sp^ ■ p + 3Jv] . (41) 

Diagonalization of the boson Hamiltonian ([|l]) is straightforward and represents an alterna- 
tive method for exactly solving the SO (8) model. 

D. Boson mean-field and fermion-pair approximations 

It is of interest to consider approximate solutions for the ground state of this model 
as well. Here, too, the natural approximations to look at first are those based either on 
mean-field boson methods or on the analogous fermion-pair approximations. 

We will not discuss in detail either the methods or the conclusions that emerge from 
these approximate treatments, since they parallel so closely the discussion for the SO (5) 
model. Rather, we will just note some of the differences that show up relative to the SO (5) 
analysis, as derived from the earlier more extensive treatment of the SO (8) model in Ref. 

One of the more interesting new features that emerges is a third type of mean-field 
solution. The solutions A and B, discussed in Sections |T| and |Tj, of course persist in the 
SO (8) analysis. But now a third solution C also appears, corresponding to a phase with 
nh — pp and T = np — ph pairs. 

As in the SO (5) analysis, there is a partial decoupling of the pairing phases in the mean- 
field solutions. Specifically, the T = 1 np — ph phase is absent from solutions A and C. This 
leads to significant deficiencies in the mean-field description of T z w nuclei, except when 
isoscalar pairing is dominant. 

E. Beyond fermion-pair correlations: A second boson mapping 

The boson mean-field and fermion-pair approximations touched on in the previous sub- 
section incorporate fermion-pair correlations. The fact that they are unable to describe the 
detailed properties of the model suggests the need for a more sophisticated procedure that 
incorporates further correlations. 

1. The T=0 S = case 

The only operators that enter the first boson Hamiltonian fl4"T|) are the three operators 
(|H) of the isovector boson space and the three analogous operators 
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pt.pt ; p. p ; pt . p (42) 



of the isoscalar boson space. The operators ( [42]) have commutation relations analogous 
to (p5|) . As a consequence, the Hamiltonian (fil]) exhibits an 0(2,1)®0(2,1) symmetry. 
Furthermore, the two-body terms of this Hamiltonian are precisely in the form of a boson 
pairing interaction. 

As in our treatment of the SO (5) model, we choose to include the effects of boson pairing 
(or equivalently four-fermion correlations) through the use of a second boson mapping. 



Furthermore, we follow the same strategy as in Sect. [TV], first carrying out the analysis for 
T = S = systems, and then generalizing. 

When dealing with 1 = 05 = systems, we must map onto a boson space defined by 
two scalar-isoscalar bosons. One is the t boson introduced in Subsect. [TVTl which reflects 



the correlations of two s bosons. The second, which we denote by q, reflects the correlations 
of two p bosons. The Dyson realization of the 0(2,1)®0(2,1) algebra contains two parts. 
The set of operators (|24"D of the isovector boson space map according to (|26|). The set of 
operators in the isoscalar space map according to the analogous relations 

p\ . p\ _> Qgt _|_ 4gtgtg ^ 

P-P ~> 9 , 

p ] • p -> -2q ] q . (43) 

Applying the second boson mapping to the Hamiltonian ([IT]) leads to a U(1)®U(1) (or 
SU(2)) Hamiltonian 

#bb = -g{i + x) [(n + - - N)th + thht - - q H - q^qt] 

5 3 
-g(l -x)[(n + -- -V)? 9 + ? ? ?? - ? - tttft ? ] 

+^ ph [6(g t t + t f g) + A(q ] q ] qt + tHHq) + 8tHq^q + 3Af] , (44) 
where now J\ is given by 

Af = 2{th + q^q) . 



In general, the problem of the Hamiltonian ( f4"4j ) must be solved by numerical matrix 
diagonalization. In the case of the dynamical symmetry limits, however, we can obtain 
simple analytic solutions. 

In the SO T (5) limit, the Hamiltonian matrix in the \t^~ nq q nq )) basis has all zeros below 
the diagonal. The same considerations as in Subsect. |IVB 2] can therefore be applied. The 
diagonal elements of the Hamiltonian matrix give directly the eigenvalues and reproduce the 
exact result @ for the T = energy when an identification n q = |t> s is made. The same is 
true for the SO s (5) limit, with the roles of the isospin and spin and also the t- and g-bosons 
interchanged. 

To obtain the solution in the SU(4) limit, we first impose the transformation 

r f = - g f ) , 



w 



V2' 

i - —(!> 4 J 



V2 



(tt + gt) , (45) 
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and rewrite the Hamiltonian (143) as 



#bb = -</[-(« + \ ~ + ~(rV - w^w) 

_l_ r t r t rr _)_ 2rVu^u> + w'w'rr] 
+4^p h [3w''w; + w'w'ww — w'w^rr] . (46) 

Working with this Hamiltonian in the basis \r^~ nw w nw )), we recover the exact SU(4) 
eigenvalues (|39"D , when an identification n w = ^\ is adopted. 



2. General case 

The procedure just outlined for carrying out a second boson mapping to the t — a space 
enables investigations of M = even T = S = states only. For a more complete treatment, 
we must appropriately extend the second boson mapping, in much the same way as we did 
in Subsection [IV B 2[ Namely, we must include, in addition to the sets of operators 
(p9f) and ([42]) , the additional one-boson creation and annihilation operators 



p ] , p . (47) 
The full second boson mapping is now given by (|3lD and 

p^ ■ p — > —2q*q + 7T ■ tt , 

pi, -»• 7f*. , (48) 

where n denotes the additional vector-isoscalar boson needed to complete the ideal boson 
space. 

This ideal boson space - formed from the t, a, q, and 7r bosons - contains an unphysical 
subspace. Since the sectors t-a and q-n are separated and since we have already determined 
the physical states (|33| ) in the t-a sector, we find that the physical basis of the full problem 
is of the form 

^i {M -T-S)-n qq n q . a T T . ^ S + T < 2VL - N . (49) 

The second boson image of the SO (8) Hamiltonian is straightforwardly obtained and in 
general can be solved by matrix diagonalization. We will discuss here only the dynamical 
symmetry limits, where analytical solutions can be obtained as in Subsection [IV B 2| . 

In the SO T (5) limit, for example, the Hamiltonian takes the form 

— cr ■ a) 

+3tH + 2thht - 2tha ] ■ a + ^ ■ ah 

I 

-3qH - 2gYgt + 2g t tvr t • tt - -n ] ■ nh] , (50) 
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with 

Af = 2{tH + g f g) - a ] ■ a - 7r f • vf . 

In the basis the exact energies (|9|) are obtained when the identification 2n q + S — \v s 
is made. The left vectors of the basis ([EJ) give the left physical eigenvectors in this limit. 

The same arguments can be applied to the SO s (5) limit. The only difference is that we 
must interchange the role of the isospin with that of the spin, and likewise the role of the t- 
and a- bosons with that of the q- and 7r-bosons. 

In the SU(4) limit, we again apply the transformation fl4"5|), after which the second boson 
Hamiltonian takes the form 

H BB = -g[-(fL + l- J\f)J\f + 3r f r + rVrr 

+2r*rw'w + w^w^rr — verier ■ a + ix^ ■ n ) 

\[2 

—nrricr ■ a — ix^ ■ 7T ) H — ~^~( a ' °^ ~ ^ ' ^ ) r ) 

+g p h{6w^w — 6r'r + Aw^w^ww — Aw^w^rr 

+3J\ — Aw^w(cx^ ■ a + 7r^ • 7f) + Aw^r(a^ ■ a — tt^ ■ tt ) 

H — -/={ a ^ ■ ^ + ^ ■ 7T^)w l={ a ^ ■ ^ — ^ ■ + 2a^ ■ an^ ■ n] . (51) 

v2 v2 

Again, the SU(4) energies (^) are easily obtained by working in the basis 

^{N-T-S)-n Ww n w . a T T . ? (52) 

with 2n w + S + T = X. The left vectors of the basis (^) give the left physical eigenvectors 
in the SU(4) limit. 



F. Wave function in terms of fermion-pair operators 

As we have just seen, the left physical eigenvectors of the second boson- mapped space 
take a very simple form in the three dynamical symmetry limits. As in Subsection [IV Q we 
can use this fact to find the corresponding wave functions in terms of fermion-pair operators. 
To do this, we must use the chains (133) as well as the analogous P-sector chains 



P ■ P — > p ■ p—> q , 

Pn -> Pn -»• • (53) 



Proceeding as in Subsect. [IV C| , we find that in the SO T (5) limit the wave functions are 
given by 



\NTT z SS z v s ) oc V Tz VsXS ] ■ St)K^-T-^)(pt . P \f 2 i\ 



Vs-S) 



x(\sl + \SU + -^Slr^P} + \pU + -^H) S \0) , (54) 
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where Vs z is the projection operator that picks states with spin-projection S z . 

An analogous expression holds in the SO s (5) limit with the roles of isospin and spin and 
also of and Pt interchanged. 

In the SU(4) limit, we must again take into account the transformation (f4~5|) . This then 
leads to 



\NTT Z SS Z \) oc V Tz V Sz (S ] ■ S ] - P f ■ P ] )^ {M - X \S^ ■ S ] + P f • P f )5 



(X-S-T) 



x(\st + \SU + ±= 2 Sl) T {\pi + \PU + ±Pt) 8 \0) ■ (55) 

In the case of the ground-state solutions, the wave functions associated with the three 
dynamical symmetry limits can be unified. For even-even N = Z nuclei, for example, all 
three can be written as 

|ATT = 0T, = 05 = 0) oc (a& - & - VP* - Pi)" /2 \0) . (56) 

In the SO T (5) limit, a = 1 and (3 = 0. In the SO s (5) limit, a = and (3 = 1. And, in the 
SU(4) limit, a = /3= l/y/2. 

Similarly, we can write the ground-state wave function for even-even N > Z nuclei by 
adding the appropriate number of isovector nn pairs to the N = Z solution (|56D , viz: 



\MT = T Z T Z S = Q) 

oc (orSt • & - /3Pt . pt)(^-r.)/2^r. | 0) (57) 



Finally, for odd-odd nuclei, we must add to (|57| ) either an isoscalar or an isovector np 
pair, depending on the symmetry limit. In particular, in the SO T (5) and SU(4) limits, the 
ground state is given by 

|A^T = T Z + 1T Z ^ = 0) 

oc (as* • 5+ - /3pt . pt)(Ar-T,-i)/ 2jS tT^t| 0) (5g) 

Analogously, in the SO s (5) and SU(4) limits, it is given byEH 

\MT = T z T z S = l) 

oc (aS-t ■ St - /?Pt ■ pt)(^-r,-i)/2 5 ^.pt| ) . (59) 

In the SOf8) model, the quartet creation operator that represents an alpha cluster takes 
Fusl 



the form 11 



14 Here, the normalization a 2 + f3 2 = 1 is used. 

15 In the SU(4) limit, there is a degeneracy in the ground states of odd-odd nuclei, explaining why 
we give two distinct ground-state solutions. 

16 There is a sign mistake in Ref. [|ll]] 
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• S* - • Pi) . 



2^/3ft(fi + 2 

We see, therefore, that in the SU(4) limit the ground state involves as many alpha-correlated 
structures as possible. 
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FIG. 3. The variational parameter a 2 in the approximate SO(8) ground-state wave functions 
(54-57) as a function of the Hamiltonian parameter x that controls the balance between isoscalar 
and isovector pairing. The behavior is shown for a) various shell degeneracies VL with fixed N = 6 
and T z = 0, b) various iV pa i r = M with O, = 10 and T z = 0, c) various T z with fixed 17 = 10 and 
J\f = T z + 6, and d) various values of the particle-hole strength g p ^ in the T = 1,S = 1 channel 
with fixed M = 6, T z = 0, and Q = 10. 

Outside the dynamical symmetry limits, the above forms for the ground-state wave 
function are not exact. Nevertheless, when we consider them as variational wave functions 
(with a and (3 as variational parameters), we find that at the minimum they have almost 
perfect overlap with the exact wave functions. Fig. |3] illustrates the behavior of a 2 as a 
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function of the parameter x of the SO (8) Hamiltonian for a few representative cases.0 

We first discuss the results in the limit g p h = 0, as illustrated in Fig. ||a-c. In this limit, 
we are in the pure isovector phase at x = 1 and in the pure isoscalar phase at x = — 1. 
Furthermore, the transition between the two phases becomes sharper as Q increases and/or 
the number of pairs decreases. In passing from the one phase to the other, we of course 
pass through the SU(4) phase. Note, however, that this is done smoothly. There is no 
plateau at the SU(4) phase, suggesting no special stability associated with the maximally 
alpha-correlated state. 

When a particle- hole force is turned on [see Fig. ^]d], some interesting changes show up. 
In particular, the pure isoscalar and isovector phases are never realized. And, furthermore, 
the transition through the SU(4) phase is less abrupt. This latter remark may have relevance 
to real nuclear systems, which are most likely near to - but not precisely at - the SU(4) 
limit. Our results suggest that even away from the SU(4) limit, there may be significant 
alpha-particle transfer strength, as long as there is a sufficiently strong particle-hole force 
present. 

Finally, we close this section by reiterating the key conclusion of our analysis of the 
SO (8) model. Namely, in the SO (8) model - a prototypical model that involves np pairing 
correlations on the same footing as nn and pp pairing correlations - the ground-state wave 
functions exhibit essential four-particle correlations. The fermion SO (8) ground state is 
constructed in such a way that the maximal possible number of fermions form correlated 
four-particle S = T = structures and the rest form like-particle pairs and/or an np pair 
and provide the isospin and spin of the system. 

VI. DISCUSSION 

We have studied neutron-proton pairing correlations within the framework of two simple, 
solvable, and somewhat realistic nuclear models. In our analysis, we made extensive use of 
boson-mapping techniques. Since the bosons introduced by the (first) mapping represent 
correlated fermion pairs, such an approach provides a quite convenient framework in which 
to investigate pairing features in fermion systems. 

With that in mind, we first studied the isovector-pairing SO (5) model, both within the 
framework of boson mean-field methods and also through the use of genuine fermion-pair 
approximations. The latter represent generalizations of the standard procedures to treat 
the pairing between like nucleons. We find that these approaches fail to describe details of 
neutron-proton pairing. In particular, the BCS solution does not allow for the coexistence 
of like-particle and neutron-proton pairs. As a result, two-nucleon transfer strengths are not 
given correctly, even though the approximate energies are quite close to the exact values. 
Tiny details of neutron-proton pairing correlations, as reflected in double binding energy 
differences, are also not well reproduced. Number and T z projection, by themselves, are 
unable to improve the situation. 



17 In the odd-odd case, the form (pq) applies for x > and the form (59) for x < 0. 
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A nice feature of the boson mapping procedure is that it provides detailed information 
on the boson-mapped Hamiltonian, information that can yield useful clues as to how to 
introduce the necessary additional correlations. The boson Hamiltonian that emerged from 
a mapping of the SO (5) fermion Hamiltonian contained an attractive pairing force between 
bosons, suggesting that correlations between pairs of bosons could be the key missing in- 
gredient to an improved description of the system. After first exploring the possible use of 
boson BCS approximation to accomodate these additional correlations, we then turned to 
an alternative treatment based on a second boson mapping. A simple description of the 
neutron-proton pairing problem was achieved in terms of the second bosons, confirming the 
importance of boson-boson (or equivalently pair-pair) correlations. Inverting the method, 
we were also able to obtain in closed form the fermion wave functions of the SO (5) model 
expressed in terms of the fundamental fermion-pair operators. The wave functions that 
emerged clearly showed four-nucleon correlations. 

We then carried out a similar study of a somewhat richer model involving neutron-proton 
pair correlations based on the algebra SO (8). We focussed our analysis on the dynamical 
symmetry limits of this model, where the group structure could be used to obtain simple 
closed expressions. Here too we found that a description limited to fermion-pair (or boson 
mean-field) correlations was insufficient and that four-particle correlations were needed to 
accurately reproduce the exact results. 

Thus, from these model calculations, we conclude that correlations involving pairs of 
fermion pairs, or alternatively quartets of fermions, are important in the regime of neutron- 
proton pairing. It is not sufficient only to pair nucleons. Whenever possible, two nucleon 
pairs will couple together to form a T = S = (alpha-particle-like) structure. 

Qualitatively, this conclusion can be understood as follows. The smallest "cluster" that 
can simultaneously accomodate two-neutron pairing correlations, two-proton pairing cor- 
relations and neutron-proton pairing correlations is one that involves four nucleons - two 
neutrons and two protons. Of course, when there is an excess of particles of a given type, 
not all particles can form these maximally-correlated alpha-like structures. Instead, they 
remain in like-particle and/or neutron-proton pairs, appended to the alpha-like condensate. 

The above conclusion is not limited, however, to cases in which all pairing modes con- 
tribute significantly. In the SO s (5) limit of the SO (8) model, for example, which only involves 
isoscalar np pairing, the ground-state wave function involves a condensate of P^ ■ P^ pairs 
and thus contains four-particle correlations. There, however, the fact that the ground state 
involves such a quartet structure is a direct reflection of angular momentum restoration. 

It is important to note, however, that the four-particle correlated structures that emerge 
are not exact alpha particles. Nor are they necessarily the most alpha-like structures avail- 
able within the model. In the SO (8) model, for example, the ground state is dominated by 
the maximal alpha structure in the SU(4) dynamical symmetry limit only. 

A challenging problem that still remains is: What is a good way to incorporate such 
four-fermion correlations into nuclear many-body approximation schemes? As noted above, 
boson BCS following a boson mapping does not seem to be an acceptable procedure for in- 
corporating four-nucleon correlations in such models; fluctuations in the particle number are 
simply too large. Another possible approach would be to include the additional correlations 
through projection. Isospin projection, in particular, would seem to be important. In the 
simple SO (5) model, for example, it leads to the exact solution. It still remains, however, 
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to develop an appropriate isospin projection technique and to apply it to realistic nuclear 
structure problems. 

We have put forth two ideas that could perhaps provide the basis for an improved the- 
ory with four-particle correlations. One possibility would be to start from a trial ground- 
state- wave function in the form of Eqs. (|56|)-(|59|). The S' and are now collective pair 
creation operators, whose structure ideally should be determined variationally (or, less ide- 
ally, from an analysis of simple two-body systems). Such an approach is a generalization 
of the generalized-seniority scheme for like nucleons, which is known to be connected to 
number-projected BCS theory. 

A second possibility worth further investigation is the use of iterative boson mappings. 
We have seen that two boson mappings, coupled with a Hartree-Bose treatment in the 
second boson space, is a way to build a number-projected theory involving four-nucleon 
correlations. In the current studies, where spurious states could be readily identified, such 
an approach proved extremely useful. Whether it will continue to prove useful in more 
realistic applications, however, where spurious states cannot be so readily separated, still 
remains to be investigated. 

Other truncation schemes in the shell model should also be studied. Our analysis of the 
SO (8) model suggests that the usual truncation schemes built in terms of separate dynamical 
symmetry limits for neutrons and protons (the SU(2) seniority limit, for example) should 
not be of much use in the general neutron-proton pairing problem. A truncation in the 
SO (5) or SO (8) seniority quantum numbers, which would then include the necessary isospin 
correlations, could perhaps be useful. 
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